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Abstract 

A systematic method is proposed for analyzing the renormalization scheme uncertainties in the next- 
next-to-leading order QCD predicitions, based on a condition which eliminates schemes that give rise 
to large cancellations in the expression for the characteristic scheme invariant combination of the 
expansion coefficients. Using this method it is shown that the QCD corrections to the tau lepton 
decay are rather stable with respect to change of the scheme, provided that an improved formula 
is used, which involves numerical evaluation of the contour integral in the complex energy plane 
with the Adler function under the integral. Optimized predictions for the tau decay corrections 
are given. It is shown that also in the case of the of QCD corrections to e + e~ annihilation into 
hadrons the conventional expansion has sizable scheme dependence, even at large energies. However, 
a considerable improvement is obtained when the QCD corrections are expressed as a contour integral, 
with the Adler function under the integral, resumming in this way the large tt 2 contributions. In the 
case of the corrections to the Bjorken sum rule for polarized structure functions it is found that for 
rif = 4 they are insensitive to change of the scheme. However, the nj = 3 expression is found to be 
strongly scheme dependent at lower energies. 
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Perturbative predictions in finite order of perturbation expansion depend on the choice of 
the renormalization scheme (RS). In the next-to- next-to-leading order (NNLO) approximation, 
when mass effects are neglected, the freedom of choice of RS may be characterized by two 
independent continuous parameters. The differences between predictions in various schemes 
are formally of higher order in the coupling constant, but numerically they are significant for 
phenomenology 

There are two things that one should do with the RS dependence of perturbative approxi- 
mants. First, using various heuristic arguments of physical or technical character, one should 
make a careful choice of the scheme, obtaining an "optimized" perturbative prediction. Several 
propositions have been discussed in the literature ]IJ , including a very interesting prescription 
based on the Principle of Minimal Sensitivity (PMS) 0. Secondly, one should investigate how 
strongly predictions change when one moves away from the preferred RS. By calculating the 
variation in the predictions over a set of a priori acceptable schemes one obtains an estimate 
of the reliability of the "optimized" prediction. A systematic method for obtaining such esti- 
mates has been recently presented in ||. This method is based on a condition which eliminates 
from the analysis the schemes that give rise to unnaturally large expansion coefficients in the 
expansion for the physical quantity or the beta-function. 

Let us consider a NNLO expression for a physical quantity 5, depending on one energy 
variable P, in the massless approximation: 

S(P 2 ) = a (P 2 )[l + na(P 2 ) + r 2 a 2 (P 2 )], (1) 

where a(/i 2 ) = g 2 (p 2 )/(4ir 2 ) denotes the coupling constant that satisfies the NNLO renormal- 
ization group equation: 

p— = -bo 2 (1 + ci<2 + c 2 a 2 ), (2) 
dp 

The expansion coefficients r« and c 2 depend on the choice of the renormalization scheme - 
the relevant formulas have been collected for example in ■ - but there exists a combination 
of these coefficients which is independent of the scheme: 

p 2 = c 2 + r 2 -ciri -r\. (3) 

As was discussed in [§], U, this combination provides a natural RS independent characterization 
of the magnitude of the NNLO correction. The RS invariant p 2 may be used to distinguish 
between "good" and "bad" schemes. Indeed, one usually identifies unnatural schemes with 
ones that introduce large expansion coefficients in an artificial way. However, the combination 
(|3]) must stay the same even for very bad schemes, which may be achieved only by the presence 
of large cancellations between various terms in the expression for p 2 [[|. It is then useful to 
introduce a function: 

o-2(ri,r 2 ,c 2 ) = |c 2 | + |r 2 | + ci|ri| + r\. (4) 

This function measures the degree of cancellation in p 2 , and gives a quantitative meaning to 
the notion of "naturalness" of any chosen RS. It is clear that a large cancellation in p 2 would 
imply that the value of a 2 would be large compared to |p 2 |. (Obviously a 2 > |p 2 |.) If we have 
any preference for using some renormalization scheme, we should also include in the analysis 
predictions obtained in schemes which have the same, or smaller, value of a 2 . For example, 
for the PMS scheme we have: 

r PMS = 0(o PMS ), c p 2 MS = 1.5p 2 + 0(o PMS ), (5) 
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which implies G2{PMS) ~ 2 [ P2 1 - More generally, we may write the condition on the acceptable 
schemes in the form: 

0-2(ri,r 2 ,c 2 ) < I \P2\- (6) 

The constant I in the condition (|6]) controls the degree of cancellations that we want to allow in 
the expression for p 2 . The value Z = 2 in (Q) is the minimal value of / for which the PMS scheme 
falls into the "allowed" region of scheme parameters. The proposition of @, H is to calculate 
the variation of the predictions for 5 over the set of schemes satisfying the condition (|6|), and 
use this variation as a quantitative estimate of reliability of the perturbative predictions. As 
was discussed in Q, it is convenient to use r\ and c 2 to parametrize the freedom of choice of 
the RS in the NNLO approximants. 

It should be stressed that no claim is made that in this way the actual theoretical error of 
the calculation is obtained, i. e. there is of course no theorem that guarantees that the difference 
between the NNLO prediction and the true result would lie within the obtained "error bound." 
However, the region of scheme parameters satisfying the condition (|6]) with / = 2 appears to be 
the minimal set of schemes that has to be taken into account — consequently, if strong scheme 
dependence is obtained for these schemes, it is an unambiguos sign that the perturbation 
series is not reliable. Also, evaluating by our method the variation in predictions for several 
quantities with the same value of I in the condition ([]) we obtain a very good estimate of 
relative reliability of the predictions. The obtained estimates may then be used for example 
to assign different weights to different observables in the global fit of Ajjg. 

Let us begin with the QCD correction 5 T to the tau lepton decay into hadrons. The QCD 
prediction for 5 T may be expressed as a contour integral in the complex energy plane ||: 

5 T = ±- r d6 (l + 2e w - 2e m - e m ) W{-°)\*=- m *eA , (7) 
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where 5n( — c") is the so called Adler function, which has the form ([!]). The contour integral 
has been initially evaluated by expanding a(— a) in terms of a(m^). This led to the frequently 
used expansion for S T in terms of a(m^), which has the form (|l|), with p\~ = —5.476 |fj, |7|. The 
region of scheme parameters satisfying the condition (^) with I = 2 is approximately given 
by ri G (—1.54,2.11) and c 2 G (—8.21,2.74). Applying to this expansion the method outlined 
above we find very strong RS dependence. This confirms observations made in M. However, 
as was pointed out in [§, |9| , one may evaluate 5 T in an improved way, using under the contour 
integral the renormalization group improved expression for 5n, and calculating the contour 
integral numerically. In this way one resumms large n 2 corrections which would otherwise 
appear in higher orders. Let us note that for 5u one has p^ = 5.238, which has the same 
magnitude, but an opposite sign compared to value obtained for the "naive" expansion. A 
complete analysis of the RS dependence of the improved expression for 5 T has been described 
in |TIJ. (0 contains some discussion of scale dependence of the improved expression.) The 
QCD predictions obtained in the improved evaluation appear to be quite stable with respect 
to change of RS, despite the low energy scale of the process. The problem of finding the PMS 
predictions for the improved expression has been considered. This is nontrivial since for the 
improved approximant we cannot use the set of algebraic PMS equations given in 0. It was 
found that the location of the critical point closest to the I = 2 region of allowed scheme 
parameters is well approximated by r\ = and c 2 = 1.5 p^, for most values of m T /A.jjs. 
Therefore predictions in this scheme have been taken as preferred predictions for the phe- 
nomenological analysis. Assuming (R T ) e xp — 3.591 ±0.036 as an averaged experimental value 
(see discussion of the experimental results in the Appendix of |H|) it was found that = 
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376(opt)ti4(th,l=2) ± 29(exp) MeV. For the / = 3 region the variation is iff MeV. This corre- 
sponds to: a^{m 2 T ) = 0.332±g;gg|(th,l=2) ± 0.015(exp). For the I = 3 region we obtain vari- 
ation of toilt- Extrapolating to m\ we find af^(m|) = 0.1 190±g;S(th,l=2) ± 0.0017(exp). 
(For the / = 3 region the variation is taooi3-) ^ e see ^ na * ^ ne ^ dependence ambiguities 
in the determination of aJf (m 2 ) or a^f (m 2 z ) from R T are not very big — this is the result 
of stabilization of the predictions when improved evaluation procedure is used. It should be 
noted however that the RS dependence ambiguities are still comparable in magnitude to the 
uncertainties related to the present experimental accuracy of R T . 

In the case of the QCD correction 5 e + e - to the the e + e~ annihilation into hadrons, for 
rif = 5, we have p| +e = —15.055 |J |7j (this does not include a very small singlet correction, 
which is added separately). The region of scheme parameters satisfying the condition @ 
with / = 2 is approximately given by r\ G (—4.8,4.2) and C2 G (—22.6, 7.5), and for I = 3 it is 
approximately r\ G (—5.5, 4.9) and C2 G (—30.2, 15.1). Changing the scheme parameters in the 
I = 2 region we find for s/ Aj^ = 75 2 the variation in the predictions from 0.0491 to 0.0537, and 
for the I = 3 region from 0.0454 to 0.0539. We see that even though the considered energy is 
high and the perturbation series should be reliable, the scheme dependence of the predictions is 
surprisingly large. This is a consequence of the fact that 5 e + e - has very large NNLO correction, 
which is reflected by the magnitude of the RS invariant. However, closer examination shows 
that major part of the NNLO correction comes from the term proportional to it 2 , which 
appears in the process of analytic continuation from spacelike to timelike momenta. The large 
contributions from the 7r 2 terms may be avoided — similarly as in the case of the tau decay - 
by expressing 5 e + e - as a contour integral in the complex energy plane, with the renormalization 



group improved expression for the Adler function under the integral [TT 

8 P + e - (s 



f d0<f n (-*)U-.e«- (8) 

lH J-TT 



For the Adler function for nf = 5 we have = —2.969. The improved expression for 5 e + e - 
appears to be much more stable with respect to change of RS as compared to the conventional 



expression [12|. Also, the I = 2, 3 regions of the scheme parameters are much smaller than 



in the case of conventional expansion, because the RS invariant is much smaller. In Fig.l 
we show the contour plot of S e + e -, as a function of scheme parameters, for s/A 2 ^ = 75 2 . 
We see that the variation of predictions over the I = 2 region of scheme parameters is in 
fact negligible from the point of view of phenomenological applications. For n/ = 3,4 the 
same effect of reduced RS dependence in the improved formula is found. Incidentally, for 
rif = 5, the optimized prediction obtained from the contour integral is close to the prediction 
obtained with the conventional expansion in the MS scheme and the PMS prediction in the 
conventional expansion. This is not the case for other number of flavors. 

In the case of the QCD correction to the Bjorken sum rule for the polarized structure 



functions |13| it was found that the expression for rif = 4, with p2 = 1.330 is very stable 
with respect to the change of the RS. The situation with rif — 3 predictions is quite different. 
For rif = 3 we have p2 = 5.476, so that the / = 2 region of scheme parameters extends 
approximately for r\ G (—1.65, 1.0) and Ci G (—2.8,8.3). Changing the scheme parameters in 
this region we find a considerable variation of the predictions for Q 2 /A 2 - T? below approximately 
6 , as is shown in Fig.2. This indicates that perturbative predictions in this region are not 
reliable, even though we may still obtain PMS predictions. Unfortunately, in the case of the 
Bjorken sum rule we cannot expect that the RS dependence would be reduced in a way similar 
to tau decay or e + e~ annihilation. Other ways of improving the QCD perturbation expansion 
in this case have to be investigated. It should be noted that the QCD correction to the Gross- 
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Llewellyn-Smith sum rule has almost identical expansion as in the case of the Bjorken sum 
rule, so that our remarks apply also to that case. 
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Figure Captions 

Fig.l. Contour plot of S e+e for s/Aj^g = 75 2 , as a function of the scheme parameters, obtained 
by evaluating the contour integral numerically. The regions of allowed scheme parameters for 
/ = 2, 3, 6, 10 are also indicated. 

Fig.2. QCD correction to the Bjorken sum rule for rif = 3, as a function of \/Q^/Ajjg, for 
different values of r\ and C2 belonging to the I = 2 allowed region: a)(-1.65,0), b)(0.5,5.27), 
c)(-1.0,5.27), d)(-1.65,5.27). The dashed curve indicates prediction in the MS scheme. 
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Fig.l 
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